LATTICE CALCULATION OF QCD 
VACUUM CORRELATION FUNCTIONS * 



M.-C. Chu(^\ J. M. Grandy(2), S. Huang^^) and J. W. Negele^^) 

Kellogg Laboratory, California Institute of Technology, 106-38, Pasadena, Califor- 
nia 91125 U.S.A. 

(2) 

^ ' Center for Theoretical Physics, Laboratory for Nuclear Science and Department of 
Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
U.S.A. 

^■^•^ Department of Physics, FM-15, University of Washington, Seattle, Washington 
98195 U.S.A. 



Submitted to: Physical Review Letters 



CTP#2113 Typeset in by Roger L. Gilson June 1992 

* This work is supported in part by funds provided by the U. S. Department of Energy 
(D.O.E.) under contracts #DE-AC02-76ER03069 and #DE-FG06-88ER40427, and the Na- 
tional Science Foundation under grant ^PHY 88-17296. 







ABSTRACT 

The first exploratory calculations of QCD vacuum correlation functions on a lattice are 
reported. Qualitative agreement with phenomenological results is obtained in channels for 
which experimental data are available, and these correlation functions are shown to be useful 
in exploring approximations based on sum rules and interacting instantons. 



PACS numbers: 12.38.Gc 
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Vacuum correlation functions for space-like separated hadron currents provide a useful 
tool for exploring the structure of the QCD vacuum which supplements our knowledge of 
hadronic ground states. These correlation functions have been determined phenomenologically 
in some channels by using dispersion relations to analyze hadron production and r-decay 
experiments/ and are also amenable to calculations using lattice field theory, QCD sum 
rules, and the interacting instanton approximation. Hence, in this work, we report the first 
exploratory lattice QCD calculation of meson and baryon correlation functions to demonstrate 
reasonable agreement with experimentally measured results where available and to show the 
potential more extensive calculations have to explore, test, and refine approximations. 

The lattice calculations were performed on a 16^ x 24 lattice in the quenched approxi- 
mation at inverse coupling ^jg^ = 5.7, corresponding to a physical lattice spacing defined by 
the proton mass of approximately 0.168 fm. This inverse coupling constant is large enough 
to give a semi-quantitative approximation and has the significant advantage that light quark 
propagators are available for point sources. Hence, propagators with point sources calculated 
by Soni et al? for 16 configurations with five values of the quark mass l/2n — l/2nc ranging 
from 351 MeV to 25 MeV were used to calculate the meson and baryon correlation functions 
listed in Table I. In order to calculate the ratio of interacting to free correlation functions, 
corresponding free lattice propagators at a very small quark mass, ma = 0.05 were calculated 
exactly on a (48)^ lattice, where the lattice volume was chosen large enough to eliminate finite 
volume effects for spatial separation less than 2 fm. Because the propagators had hard-wall 
boundary conditions in the time direction, all correlation functions were calculated on the 
central time slice containing the source. 
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Several lattice artifacts had to be eliminated to obtain a physical approximation to the 
continuum correlation functions. Significant anisotropy is introduced in the rotationally in- 
variant continuum correlation functions by the Cartesian lattice. Although at very short 
physical distances the anisotropics in the numerator and denominator cancel by asymptotic 
freedom and at a sufficiently large number of lattice spacings the granularity of the lattice 
becomes negligible, for the present lattice spacing there remains an intermediate region in 
which some correlation function ratios display substantial anisotropy. As discussed and jus- 
tified in detail in a subsequent article,^ since the lattice points far away from the Cartesian 
directions agree with the continuum result in the non-interacting case and should be the most 
reliable in the interacting theory, the lattice results in this work are taken from sites r within 
a small angular cone surrounding the diagonal directions d = (n, n, n) such that d • r > 0.9. 
Also, although in principle one would like to normalize the ratios at an infinitesimally small 
separation, we have had to normalize our results at the first non-zero diagonal separation 
(1, 1, 1), corresponding to a physical separation of ~ 0.29 fm. 

A second lattice artifact arises from the contributions of periodically repeated images of 
the physical sources on the finite lattice having periodic boundary conditions. As discussed 
in detail in Ref. [4] , it is straightforward to perform a self-consistent subtraction of the image 
contributions, and we have applied and verified this image correction to all the correlation 
functions presented in this work. 

Finally, because it is impractical to calculate quark propagators at a quark mass light 
enough to produce a physical pion, it is necessary to perform a sequence of calculations at a 
series of heavier quark masses and extrapolate to the pion limit. Although we do not know 
the proper functional form for the extrapolation in the chiral limit, the extrapolation is in fact 
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innocuous because we only have to extrapolate data at ruq = 351, 199, 110, 67, and 25 MeV 
down to 8 MeV. A worst case example will be shown below in Fig. 1. Thus, we believe the 
technical aspects of relating the present lattice calculations to the physical continuum theory 
are under reasonable control, and that our results provide a meaningful first comparison of 
quenched QCD with phenomenological and model results. 

The correlators we considered are listed in Table I. For convenience, we have contracted 
the vector indices in the vector and axial channels. In addition, we inserted a factor of Xf^'y'^ in 
the baryonic correlators before taking the Dirac trace in order to project out the component 
which is stable in the massless quark limit. For a local field theory such as QCD, a two- 
point function is uniquely characterized by its absorptive part in momentum space (up to a 
polynomial) through the dispersion relation 



where the upper, middle and lower terms in brackets refer to scalar and pseudoscalar mesons, 
vector and axial vector mesons and baryon channels respectively and terms which vanish 
for the correlation functions in Table I have been omitted. Phenomenologically, we expect 
f{s) has two major contributions, a resonance piece and a continuum piece, parameterized as 
f{s) = A^(5(s — M^) + fp{s)9{s — So), where M is the resonance mass, A is the coupling of the 
current to the resonance state and fp{s) is the perturbative contribution to the correlator. 
Since QCD is asymptotically free, the /p(s)'s are well-approximated by the corresponding 
free correlators with massless quarks, also listed in Table I. This type of parameterization is 
widely used in QCD sum rule calculations. 




(1) 
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An inverse Fourier transform of Eq. (1) defines the foUowing phenomenofogical correlators 
in coordinate space as a function of M, A and sq, where the polynomial P{q^) only contributes 
at the point x = and can be ignored for finite x: 

3M^x-^ I ^/ 1 \ J 1 1 ^ (2) 



Here, m refers to the upper terms in brackets for scalar and pseudoscalar mesons and the 
middle terms in brackets for vector and axial vector mesons, and b refers to the lower terms 
in brackets for baryon channels. In practice, we always normalize R{x) by the corresponding 
free correlator with massless quarks R'o'ix) oc x' and R^{x) 

In Fig. 1 we show the complete lattice data and the three-parameter fit in the pseudoscalar 
meson (pion) channel for each of five values of the bare quark mass. The result at the physical 
pion mass is obtained by binning the lattice data in two-lattice-unit bins and extrapolating 
the binned data with the results shown by the solid circles. The striking result in this channel 
is the extremely rapid rise in the correlation function ratio, necessitating a log plot, which 
arises from the strong attraction and corresponding light pion mass. The lattice result agrees 
qualitatively with Shuryak's phenomenological estimate,^ denoted by the dot-dashed line, 
based on the value = (480 MeV)^ and the fact that the peak is proportional to Att/to^ ~ 
fn/'nT'q- explains the particularly large quark mass dependence in this channel. Detailed 
treatment of the extrapolation, error analysis of fitted parameters, and lattice renormalization 
corrections will be deferred to the longer paper, ^ and we only show extrapolated results and 
phenomenological fits for all other channels, where the quark mass dependence is much weaker. 

Figure 2a displays the vector meson channel (rho) result. As emphasized by Shuryak, 
the salient feature in this channel is the fact that although the free correlator falls four orders 



of magnitude between 0.3 and 1.5 fm, the ratio is nearly one over the whole range, and our 
lattice result is consistent with his phenomenological analysis of e+e~ — > even number tt's, 
denoted by the dashed curve. The result in the axial meson channel, (^i) shown in Fig. 2b 
is qualitatively similar to the phenomenological analysis of r — > Stt decay ^ denoted by the 
dashed line, although finite lattice effects render it difficult to reproduce the rising tail due to 
mixing with the pion. The result for the scalar meson channel, for which the extrapolation 
was more problematic than any other, is shown in Fig. 2c. Note that in both the axial and 
scalar channels, Eq. (2) did not produce reasonable physical parameters when fit to the data, 
so the solid curves are smooth curves to guide the eye in these cases. For comparison, the 
predictions of the interacting instanton approximation^ for mesons using a Pauli-Villars cutoff 
Apv = 130 MeV are shown by dotted lines, and in each case are in qualitative agreement with 
the lattice results. Comparison of calculations on the same lattice with "cooled" configurations 
which are in progress will be particularly instructive in this connection. 

Since there are no phenomenological data in the baryon channels, in Fig. 3 we have 
compared nucleon and delta correlators with the sum rule calculations of Belyaev and loffe^ 
(dashed lines) and Farrar et alJ (dot-dashed lines). Given the agreement with phenomenol- 
ogy in the meson channels and drastic differences in sum rule results for the delta, lattice 
correlation functions can play a useful role in testing and refining sum rule approximations. 
In summary, we believe these exploratory calculations demonstrate the feasibility and utility 
of calculating vacuum correlation functions on a lattice and motivate definitive calculation on 
larger lattices with >6. 

It is a pleasure to thank Edward Shuryak for stimulating discussions at the Aspen Center 
for Physics and elsewhere and for making his data available to us prior to publication. We 
also thank Amarjit Soni for making his point propagators available to us, and the National 
Energy Research Supercomputer Center for Providing Cray-2 computer resources. 
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FIGURE CAPTIONS 

Fig. 1: Ratio of interacting to free correlation functions measured at five quark masses (open cir- 
cles), binned data extrapolated to the physical pion mass (closed circles), three-parameter 
fits (solid lines) and the phenomenological results of Ref. [1] (dot-dashed line). 



Fig. 2: Extrapolated ratio of meson correlation functions (closed circles) and fits (solid curves) 
as in Fig. 1. Dashed lines denote phenomenological results^ and dotted lines show the 
interacting instanton approximation.^ 



Fig. 3: Extrapolated ratio of baryon correlation functions (closed circles) and fits (solid curves) 
as in Fig. 1. Dashed and dot-dashed lines denote sum rule results from Ref. [6] and [7], 
respectively. 
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Table I: Hadron currents and correlation functions. 


Channel 


Current 


Correlator 




Vector 


,/,, = Ti,'y,,d, 


(O \TJ„(x)Jn(0)\ O) 


1 

47r2 


Axial 






1 

47r2 


PspiidoscalaT 




lTJPfx)J^'fO)l 0\ 

\ V7 1 -1. I %AJ ) tJ \ / 1 / 


3s 
87r2 


Scalar 


J* = ltd 


(O TJ*(a;)P(0) O) 


3s 
87r2 


Nucleon 




^Tr ((0|TJ^(a;)J^(0)|0>x^7^) 


S2 
6477^ 


Delta 


= eabc [«"C7^w'] 


((0 |TJ^(x)J^(0)|0>x.7l 


3^2 
2567r4 



9 



